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Integration of the next generation of photonic structures with electronic and optical on-chip com- 
ponents requires the development of effective methods for confining and controlling light in subwave- 
length volumes. Several techniques enabling light coupling to sub-wavelength objects have recently 
been proposed, including grating- [13], and composite-based solutions [4 8]. However, experimental 
realization of these couplers involves complex fabrication with ^ lOnrn resolution in three dimensions. 
One promising alternative to complex coupling structures involves materials with vanishingly small 
dielectric permittivity, also known as epsilon-near-zero (ENZ) materials[9, 10]. In contrast to the 
previously referenced approaches, a single fiat layer of ENZ-material is expected to provide efiicient 
coupling between free-space radiation and sub-wavelength guiding structures. Here we report the first 
direct observation of bulk-ENZ-enhanced transmission through a subwavelength slit, accompanied by 
a theoretical study of this phenomenon. Our study opens the door to multiple practical applications 
of ENZ materials and ENZ-based photonic systems. 

It has been known that ENZ materials possess some unusual properties [9]. In particular, the electromagnetic 
field inside the ENZ structures tends to become completely homogeneous, reflecting the dramatic extension of local 
wavelength caused by a vanishingly small refractive index. More recently, it has been shown theoretically that 
lossless ENZ systems enhance coupling between two planar waveguides through an ultra-thin guiding channel. Ex- 
perimental verification of this principle, in the microwave regime, was performed in systems where the ENZ coupler 
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FIG. 1; Bulk properties of ENZ material; Panel (a) provides a schematic of the optical measurements on the bulk ENZ 

material, and the related transmission and reflection data from these measurements; symbols show experimental reflection 
(red circles) and transmission (blue diamonds); solid lines represent theoretical fits (see Methods); panel (b) shows retrieved 
permittivity spectra. 
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FIG. 2; Geometry of the system; (a) SEM image of a representative ENZ-coupled subwavelength aperture used in exper- 
iments (length of the yellow scale bar is 6/im). (b) schematic drawing of the structure geometry. The total thickness of the 
ENZ material, n-doped InAsSb, is di + dg ~ 1.15/im. The thickness of both the gold (and the slit) is dg ~ 300nm. The slit 
width (w), for the samples studied, varied from l.O/xm to 1.8/im in 0.2/xm increments. 



utilized split ring resonator structures or alternatively, is replaced with a waveguide designed to mimic the optical 
properties of an ENZ material [l2l| . While the original waveguide coupling structures discussed in the literature are 
not directly applicable to macroscopic plane- wave-to-nanoscale couplers, an alternative mechanism for ENZ-enhanced 
transmission of plane waves through a relatively thick ENZ cover, discussed using a ray optics "focusing" formalism, 
was presented in [l3|. However, to the best of our knowledge, up to now no experiments with homogenous bulk ENZ 
materials have been performed. In this work we present the first experimental study of ENZ-enhanced transmission 
through subwavelength features with a bulk ENZ material. Our material is semiconductor-based and inclusion-free, 
and operates at optical frequencies (in the mid-infrared wavelength range) , opening the door to integration of ENZ 
materials with cutting edge optoelectronic materials and devices. We also present a comprehensive analysis of this 
unique phenomenon and discuss the relationship between two ENZ-enhanced transmission geometries. 

The ENZ material used in this experiment consists of the heavily doped, narrow band-gap semiconductor, 
InAsajSbi-a;, with x — 0.89. Our InAsSb layer was grown by Molecular Beam Epitaxy (MBE) in a Riber 32 system, 
on a semi-insulating GaAs substrate. The InAsSb epilayer is 1.15/xm thick and heavily n-doped ('-^ 1 — 2 x lO^^cm"^). 
It is well known that conduction band electrons in semiconductors can behave as a free-electron plasma, whose in- 
teraction with incident radiation fields can be described by a Drude-like dispersion model. In such a system, the 
permittivity of the semiconductor can be described using: 

where w is the radial frequency of the radiation, F is the damping constant, e^o is the high-frequency permittivity 
of the semiconductor lattice, ujp is the semiconductor plasma frequency, n is the carrier concentration, and m* the 
effective mass of the free electrons in the semiconductor. By using the InAsSb material system, which can be doped 
heavily n-type, and which has a low conduction band electron effective mass, we can push the plasma frequency of 
the material well into the 6 — lO/im wavelength range. 

The bulk optical properties of our InAsSb layer have been derived from optical measurements as described in 
Methods, and demonstrate that the material exhibits ENZ behavior in the vicinity of Ao — 8/im. The resultant 
spectral dependence of the real and imaginary parts of the InAsSb permittivity are shown in Fig|TJ Values of e 
for our material are of the same order as e for noble metals (Ag, Au) or more exotic materials such as transparent 
conducting oxides [14] about their plasma frequencies in the near UV and near-IR wavelength ranges, respectively. 
Therefore, our studies provide an outlook for the performance of realistic plasmonic materials across the optical 
frequency ranges, from UV to mid-IR. In order to experimentally demonstrate ENZ-enhanced transmission through 
subwavelength optical features, subwavelength slits in an optically thick metal layer were fabricated onto the InAsSb 
epilayer, as shown in Figl2] (see Methods). At the same time, a set of control samples, with slit widths identical to 
the InAsSb epilayer samples, were patterned on a semi-insulating GaAs substrate, in order to compare the ENZ- 
enhanced transmission to structures with identical geometries, but without material losses. The GaAs substrate, 
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FIG. 3: Analysis of ENZ-enhanced transmission; The four panels represent the analysis of transmission through ENZ- 
enhanced and control samples; Panels (a,b) show transmission spectra for TM (red) and TE (magenta) polarized light, for the 
systems with (a) w = l/im and (b) w = lAfim; Panel (c) illustrate the transmission through the control structure where GaAs 
replaces ENZ material {w = Ifim); Panel (d) shows the transmitted power at the ENZ-enhanced transmission peak A ~ 8fim 
as a function of slit width, and demonstrates one of the essential characteristics of ENZ-enhanced tranmission: its extremely 
weak dependence on the slit size w. In all panels, symbols and solid lines represent experimental data, dashed lines represent 
the results from the analytical technique developed in this work; The relatively poor agreement between the analytical results 
and experimental transmission through GaAs control samples is explained by relatively large effective slit width (nw) in these 
systems. 



having relatively large refractive index (n ~ 3.3) represent the best-case scenario for traditional direct coupling of 
light into the small slit. 

The transmission spectra of the ENZ-based and the control structures are shown in Figl3lja-c). It can clearly be 
seen that the InAsSb samples provide a drastic enhancement for TM-polarized light coupling into the slit structure 
at or near the ENZ wavelength. In fact, for systems with narrower slits ('^ l/im), a higher transmission intensity 
predicted (and experimentally observed) for the ENZ samples (at Aq ^ Sfim) then their dielectric high-permittivity 
control counterparts, despite the significant losses associated with free carriers in the ENZ samples. Note that the 
ENZ material does not provide transmission enhancement for TE-polarized radiation. 

To understand the origin of the observed phenomenon, transmission through the system has been modeled with 
mode- matching software, developed according to Ref . [l5| , followed by the development of an analytical technique 
describing transmission and reflection by the system, outlined in Methods. These studies clearly demonstrate that 
transmission through the system is dominated by the fundamental (fcx — 0) mode of the planar guide formed by the 
slit. Since this mode only exists for TM-polarized light, TE waves do not exhibit ENZ-related transmission peaks. 

The amplitude of the fundamental mode, and correspondingly, the total light transfer through the system, is 
related to the interplay between two processes: (1) coupling between the incident radiation and the = mode at 
the entrance slit of the device; and (2) propagation of the mode through the slit. The flat ENZ layer at the entrance 
side of the device drastically enhances the coupling between free-space and the slit (guide). This enhancement is 
related to amplification of the amplitude of the electromagnetic field inside the ENZ material [lB| ■ In our system, this 
enhancement originates from bulk plasma resonance of the doped semiconductor. 

The secondary enhancement comes about from modal inter-coupling between the ENZ-filled guide and the bulk 
ENZ coupling layer. However, the extra in-coupling efficiency, provided by the ENZ- layer inside the waveguide is 
more than compensated for by the negative effect of the ENZ-related loss of the device as the mode propagates through 
the system, which ultimately reduces the amplitude of the fundamental mode at the exit side of the slit. Since the 
fabrication of ENZ-filled systems is significantly more complicated than the relatively straightforward fabrication of 
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FIG. 4: ENZ-assisted light funneling at different losses Panel (a) represents the comparison between transmitted power 
at the ENZ wavelength in three ENZ-assisted structures: (i) the system realized in our experiments (red), when ENZ material 
fills the substrate and in-slit regions, (ii) the system where the ENZ material is deposited as a thin substrate layer, but air fills 
both the slit region and the region above the slit (green), and finally, (iii) the design originally suggested in Ref.[13] where ENZ 
fills substrate, slit, and superstrate regions of the geometry shown in Fig[2] (black) . It is clearly seen that at relatively large 
absorption, the winning design calls for low- loss material (air) inside the guide and above it; the practical benefits of ENZ are 
limited to extending the wavelength inside the system. Lower-loss systems, however, benefit from the extra ENZ layers inside 
and behind the slit; In the limit of low loss, the ENZ light collimation, suggested in Ref.[l^ is seen in the structure with the 
ENZ superstrate. Panels (b-d) illustrate field distribution in the three systems (i-iii, respectively) at imag{e) = 0.1. 



subwavelength geometries above planar ENZ materials, the benefits of complex systems have to be weighed against 
their expected performance. Such a comparison is shown in FiglH It is seen that for realistic values of loss, transmission 
through the subwavelength channel is maximized for structures where the only ENZ component is the thin coupling 
layer, and the slit is free of ENZ material. When losses are substantially reduced, however, the positive in-coupling 
boost provided by the ENZ material inside the slit, accompanied by an ENZ out-coupling layer, add performance to 
the system. 

To conclude, we have reported the first comprehensive study of ENZ-enhanced light transmission through an 
isolated subwavelength slit. An analytical model describing the transmission has been developed, and the interplay 
between different contributions to ENZ-enhancement have been analyzed. Theoretical results are in agreement with 
experimental data. The developed formalism opens the door for ENZ-enhanced coupling to a variety of subwavelength 
plasmonic, beamsteering, nanophotonic, and nanoelectronic jj, @, 0, [l7i - t22l] systems across a broad range of optical 
frequencies. 

The authors gratefully acknowledge N. Engheta (UPenn) for useful discussion and comments. The SEM image 
in Fig. H was taken by Dave Safford at lES Technical Sales with a PhenomWorld SEM model PHENOM. This 
work was supported by the AFOSR Young Investigator Program under Grant #FA9550-10- 1-0226 and NSF Grant # 
ECCS-0724763. 



METHODS 



Permittivity of doped semiconductors The bulk optical properties of our ENZ layer can be determined from 
the fitting process outlined in Ref. |23l], specifically developed for a narrow band semiconductor (InAs) deposited 
upon GaAs. Transmission and reflection data were collected on the as-grown InAsSb/GaAs wafer, and this data 
fitted by adjusting the plasma frequency cjp and damping constant F until excellent agreement between modeled and 
experimental data was obtained. From this, the plasma wavelength of the material was determined to be 8fim 
(^ 1250cm~^) and the damping constant F — 50cm~^, in good agreement with the results of Ref. [isj. It should be 
noted that in Fig. [3] the plasma frequency does not, in fact, correspond to the minimum (maximum) of the reflected 
(transmitted) signal, as this data is a convolution of the InAsSb material properties and the Fabry-Perot cavity formed 
by the InAsSb layer on the GaAs substrate. In addition, the theoretical transmission data depicted in Fig. |3]has been 
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scaled by a loss factor exp(— 27rZo/A), with Iq ~ 1.1/zm that describes absorption from lattice-mismatch- induced defects 
in the transition layer between the GaAs substrate and the ENZ layer. Alternative models of permittivity that do 
not involve losses in this transitional layer yield qualitatively incorrect results (see supporting online documentation) 

Device Fabrication. The device was prepared using standard lithographic techniques. Five columns of slits, 
with slit widths ranging from l/im to 1.8/im in 0.2/im increments, were patterned onto the semiconductor substrate 
as a photoresist etch mask. The samples were then etched 300nm using a citric acid:hydrogen peroxide etchant. 
The photoresist etch mask is retained, and the sample coated in 295nm of Au, following a 5nm Ti adhesion layer, 
after which the slit feature is lifted off. Both InAsSb epilayer samples and undoped GaAs samples were fabricated 
simultaneously, so that transmission through ENZ and non-ENZ slit structures could be compared. 

Device Characterization Slit transmsission data was collected on a Bruker IRl infrared microscope coupled to 
a Bruker V70 Fourier Transform Infrared (FTIR) Spectrometer, with a spectral resolution of 32cm~^. The internal 
FTIR globar source acts as the incident radiation source, and a wire-grid polarizer, placed in the incident beam 
path, directly before the sample, allows us to alternate between TM (electric field perpendicular to the slit) and TE 
(electric field parallel to the slit) polarized incident light. The microscope allows seamless switching between optical 
microscopy and FTIR spectroscopy, allowing for accurate alignment of the detection optics over the subwavelength 
slit, which is clearly discernable in the visible optical microscope. Ultimately, the transmitted radiation is collected by 
a liquid nitrogen cooled HgCdTe (MCT) detector, with a 16/im cut-off wavelength. Data was collected for a spectral 
range from 700cm~^ to 4000cm~^. Figure [3] shows the TE and TM polarized transmission for Ifim and 1.4/zm slits 
fabricated on both the InAsSb and GaAs material. All transmission data was normalized to transmission through 
an unpatterned GaAs substrate to remove the effects of substrate absorption, as well as the spectral shape of the 
emission of the mid-IR internal FTIR source. A clear transmission peak is seen at ENZ for all TM polarized data on 
the InAsSb samples, with no such peak observed for the GaAs samples, in either polarization. Comparing the InAsSb 
slit transmission spectra to the unpatterned InAsSb transmission spectra shown in Fig. [l] it is clear that the ENZ 
peak is not related to any resonant transmission peaks in the bulk InAsSb material. 

Mode-matching calculations provide unprecedented insight into the optical properties of complex structures. 
In this approach, the system under investigation is separated into a set of bulk layers and a waveguide layer. The field 
inside bulk layers is represented as linear combination of plane waves (free-space modes), parametrized by x-component 
of their wavevector. The field inside the slit (waveguide) region is represented as a sum of guided (waveguide) modes 
with their field concentrated inside the slit, accompanied by a set of free-space-like bulk modes with the field extending 
throughout the metallic cladding. The spectrum of the bulk modes in the waveguide region of the structure matches 
the spectrum of the plane waves in the free-space-like regions. 

Gontinuity of the and Hy components of the electromagnetic fields is used to deduce the set of boundary 
conditions relating the amplitudes of the modes in the adjacent layers. The implementation of the mode-matching 
simulations used by our group is described in Ref. 1^ . The comparison between mode matching simulations and other 
numerical techniques for solving Maxwell equations is presented in Supporting Online Information. 

Analytical calculations of light transmission through the system. To calculate transmission through a 
subwavelength slit, we represent the incident field as a plane wave with a fixed wavevector k^; (fc^ = is used 
in the calculations), and a given amplitude Jo; and approximate the refiected field as a combination of the plane 
wave with the same wavevector and amplitude i?o, and a scattered wave with electric field represented as = 
S J e^^'^smc{wkx /2)dkx; further, the transmitted field is approximated a.s E^ — T J e^^'^smc{wkx /2)dkx (the spectral 
shape of the reflected and transmitted fields represent the field diffracted by a single slit in the perfect metal conductor 
and are consistent with results of our numerical calculations). Finally, we use the continuity of E^ and Hy along the 
slit entrance and impedance boundary conditions f23] along the metal surface to deduce the values of the parameters 
Ro, S, T, and the amplitudes of the fundamental guided mode inside the slit in terms of the parameter Iq- 



The overlap- integral formalism, outlined in Refs.[15|, [251, l26|, yields the following relationship between the parame- 
ters: 

Io + Ro + S^ = A+ + A-<!>y 
j^[Io -Ro-S- Ik{w,ei)] = ^[A+ - 

Ro ^ IoP + S ■ Is{w,ei,em,q) (2) 
'-^T = Aj^, + A- 
^T-h{w^e2) = ^[A+%~A-] 

where ei,eg,em,e2 describe permittivities of materials in front, inside, around, and behind the slit (Figl2]), ka{kx) = 
\J EaUj'^ /c? — kj^ , k^ = ka{k^), $g = exp^ikgdg), Fresnel coefficient p — ^^^""^"'^i , and the integral relationships are 
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FIG. 5: ENZ-Accuracy of numerical simulations and analytical resuls. The figure shows the agreement between the 
light transmission in the three ENZ-assisted structures described in FigU) Crosses, squares, and solid hues correspond to results 
from FEM simulations, mode-matching solutions of Maxwell equations, and our analytical formalism (Eq.©), respectively; 
Colors are identical to FigU) The results of all FEM models are scaled by the same number to match the results of mode- 
matching and analytical calculations, which are plotted with no scaling. 



given by: 

Ih{w,ea) = / TTiH ^ 

J —oo '^ay'^x } 

I {w e e q) - ^a^Pl /"°° ^gkpjkx) + epkgjk^) k^ sin{k^w/2) - q cos{kxw/2) 
fafc" + e^fcS y_oo ka{kx)kp{k,^) kl+q^ 

The parameter g ^ 1 plays the role of the fitting parameter, fine-tuning the role of the impedance boundary conditions. 
Our studies indicate that the numerical values of the results only weakly depend on the value of q (see Supporting 
Online Materials). Here we use q — 0.1. 

To account for the effects of the other interfaces (GaAs/ENZ, superstrate/air) either above or below the sub- 
wavelength feature in the systems analyzed, Eqs.(l2|) are incorporated into a conventional plane- wave transfer-matrix 
formalism [241 . 



SUPPORTING ONLINE MATERIALS 



Accuracy of Numerical Simulations The accuracy and stability of our mode-matching simulations were tested 
for convergence with respect to variations in spectrum of bulk modes, as well as the number of guided modes included 
in the simulations. Once optimal spectral representation was identified, results of mode-matching solutions of Maxwell 
equations were compared against the results of a commercial finite-element solver [23| . The agreement between the 
transmission spectra calculated using the two techniques are shown in FiglS] 

Alternative models for permittivity, and scaling of analytical data 

As explained in Methods, the Drude model tends to overestimate the transmission through the as-grown ENZ 
material, indicating the presence of additional losses in the ENZ structures. The origin of these losses can be two-fold. 
First, because the lattice constant of InAsSb does not match that of the GaAs substrate, there will be a transitional 
layer at the interface between the two materials that may result in additional transmission losses in ENZ structures. 

Alternatively, the "extra" losses could be due to impurities in the ENZ material and thus could be uniformly 
distributed through the ENZ layers, causing significant deviation of the dispersion of this layer from Drude model 
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wavelength (iim) wavelength (jim) wavelength (^m) 

FIG. 6: Alternative models of ENZ permittivity. Panel (a) shows permittivity of InAsSb, extracted under the assump- 
tion that additional loss is uniformly distributed through the structure; panel (b) shows theoretical (lines) and experimental 
(symbols) transmission and reflection spectra of a homogeneous layer of ENZ on GaAs substrate. Finally, panel (c) shows the 
comparison between ENZ-assisted transmission through w = l^m slit, measured experimentally (solid line), and calculated 
with Drude-based model (crosses) and with the permittivity model shown in panel (a) (stars). 

predictions. To assess the applicability of these two competing models, we have used a nonlinear fitting approach 
to extract the permittivity of InAsSb layer directly from tranmission/reflection measurements (not relying on Drude 
model). The results of this procedure are presented in Fig|6l[a,b). 

Furthermore, we calculated the predicted transmission spectra according to the two models of permittivity and 
compared the results to the experimental data. These calculations are shown in FigEl^c). It is clearly seen that only 
Drude-based model is consistent with experimental results. We therefore conclude that the extra loss seen in our 
system is localized in a thin transitional layer between InAsSb and GaAs materials. In accordance with this model of 
our material, the results of our analytical calculations for light transmission through the systems with ENZ substrates 
in Figs 1TJ-|4] are scaled by the loss factor exp(— 27rZo/A), with lo 2± 1.1/im. 
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